Supercurrent transferring through c-axis cuprate Josephson junctions with thick 

normal-metal-bridge 
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With simple but exactly solvable model, we investigate the supercurrent transferring through the 
c-axis cuprate superconductor-normal metal-superconductor junctions with the clean normal metal 
much thicker than its coherence length. It is shown that the supercurrent as a function of thickness 
of the normal metal decreases much slower than the exponential decaying expected by the proximity 
effect. The present result may account for the giant proximity effect observed in the c-axis cuprate 
SNS junctions. 
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In a superconductor-normal metal junction, it is con- 
sidered that the Cooper pair can penetrate into the 
normal metal within a distance of coherence length 
due to the proximity effect^ Therefore, according to 
the proximity theory, the supercurrent cannot transfer 
through a superconductor-normal metal-superconductor 
(SNS) junction when the length of the normal metal 
is much longer than £„. However, the supercurrent 
in high-temperature-superconductor (HTSC) junctions 
with very thick barrier (consisting of weakly doped non- 
superconducting cuprates) has been observed by a num- 
ber of experimentsi£i2i^£&£&2 There have been some the- 
oretical explanations based on the assumption of the exis- 
tence of superconducting puddles in the pseudogap states 
of the cuprates i 10 ' 11 But the physics of pseudogap states 
of the cuprates is not clearly understood so far. The 
problem whether the supercurrent can transfer through 
a long bridge of normal metal between two superconduc- 
tors is still an outstanding puzzle. 

The supercurrent stems from the motion of paired car- 
riers in the superconductor. It is known that the su- 
percurrent can be conducted by Andreev reflections in 
the SNS junctionsia > Miifi5 1 i6 1 i7 1 i8 1 i9 1 20 1 2i,22 1 23 The su _ 

percurrent in one superconductor, for example in the left 
one, can transfer through the SN interface by generat- 
ing the propagations of electrons and holes in the normal 
metal due to the Andreev refection. At another NS in- 
terface, the electrons and holes are converted into paired 
electrons in the right superconductor i2£ As a result, the 
paired particles are conducted from the left to right su- 
perconductors even though Cooper pairs cannot survive 
in the normal metal. In case of clean normal metal 
without large damping in particle propagations, the su- 
percurrent may transfer through the long SNS junction. 
Most of the existing calculations were based on various 
approximations including the Andreev and WKB ones 
neglecting the normal reflections ! 13 ' 14 ! 15 ' 16 However, the 
Andreev approximation is not adequate for the bound 
states of energy within the superconducting gap^ 2 - For 
long normal metal junctions, the contribution from the 
bound states to the supercurrent is significant because 
the number of the bound states is proportional to the 



length of the normal metal. Moreover, for long normal 
metal SNS junctions, one needs to carefully deal with 
the wavefunctions of the quasiparticles since the super- 
current sensitively depends on them. A reliable study 
on the supercurrent in the long SNS junctions is is still 
lacking. 

In this work, on the basis of the Andreev-reflection 
approach, we study the supercurrent in the c-axis cuprate 
SNS junctions using simple but exactly solvable model. 
We will show that the supercurrent can transfer through 
the junctions with the normal metal much thicker than 
its coherence length. We intend to provide a possible 
explanation to the relevant experiments. 

We consider a c-axis cuprate SNS junction with the 
normal metal occupying the layers from 1 — I to I — 1. A 
sketch of the junction is shown in Fig. 1. Within the ab- 
plane, the quasiparticles are described by the t — t' tight- 
binding model with t' /t = —0.3. The phase difference 4> 
between the pair potentials of the two superconductors 
drive the supercurrent. With an unitary transformation 
U — exp(i<73(/>/4), one can show that the physical quan- 
tities of the system depend only on the total phase dif- 
ference of the pair potentials. For convenience, we here 
set the phases of pair potentials as ±0/2 respectively for 
the left and right superconductors. The electron motion 
along c-axis is described by the interlayer hopping t c with 
t c /t << 1. The magnitude of t c may be the same order 
as A. The parameter of electron hopping through the SN 
interface is to- Through out this paper, we use the units 
of e = Ti — 1 with — e as the charge of an electron. 

The states of the quasiparticles are described by the 
Bogoliubov-de Gennes (BdG) equation^ 5 - Since the mo- 
mentum parallel to the interfaces is conserved during the 
motion of the quasiparticles through the junction, the 
transverse (orthogonal to the c axis) part of the wave 
function can be taken as plane waves. The problem is 
then reduced to solve the one-dimensional BdG equa- 
tion along the z-direction. For an eigenstate of trans- 
verse momentum k±, the chemical potential in the BdG 
equation is then substituted with p-(kx) = fi — e(fcj_) 
where e(k±) — — 2i(cosfci +cosfc2) — At' cos k\ cos is 
the in-plane single-particle energy with k\ and ki the two 
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FIG. 1: (Color online) Sketch of a c-axis cuprate SNS junc- 
tion. The interlayer hopping is described by t c . to is the 
hopping through the NS interface. 



components of k±. The order parameters of the super- 
conductors are then given by A(k±) exp(±icf>/2) (+ and 
- for left and right superconductors, respectively) with 
A(fcj_) = A(cos/ci — cosfc2). The BdG equation reads 



with 



Hi-; = 



^Hijtpij) = Eip(i), 



[V{ - p,(kx)]aa + A l cr+ + A*cr , for i = j 
-t c (T3, for nearest — layer hoppings 
— to 03, for interface hoppings 



(1) 



where Vi = Vq for 1 — I < i < I — 1 or otherwise, Aj = 
A(fcjJ exp(i(j)/2) for i < -I, A* = A(fcj_) exp(-i0/2) for 
i > I, and er's are the Pauli matrices. The potential 
shift Vq controls the density difference between the nor- 
mal metal and the superconductors. All the states in a 
complete basis can be divided into three types: incoming 
waves of free states from the left and right superconduc- 
tors, and the bound states mainly confined in the normal 
metal with damping tails in the two superconductors. 

The free state with an incoming wave number k + from 
the left superconductor is obtained as^ 

M) = te)( e ifc+2 '+6 e - ifc+z ') + ateu ifc 



z' = j + 1 < 

W) = lR P -iwxR„P« ' 1-I<J<1-1 



z' = j - I > 



with the boundary conditions 



*(-0-feH) = o 
^1(1-0-^/2(1-0 = 

r^ 2 (i-l)-^3G-l) = 
^2(0-^3(0 = 



where = exp(i</>/4), the wave numbers gi, 92, fc + and 
fc~ satisfy the equations £(<?i) + Vo = — £(#2) — = 
V^ 2 (fc + ) + A 2 (fc_L) ee E k and £(Ar) = H(fc + ) with 
= -2t e cosfc - /i( fe ± ), 71 = ^1/2 + e(fc+)/2^ fc , 
V = ^1/2- £(k+)/2E k , and r = t /t c . The incoming 
wave number k + is defined in the ranges — fco < k + < 
and ko < k + < tt, where the group velocity dEt/dk 
is positive, with kg = arccosf— fj,(k±)/2t c ] as the 'Fermi 
wave number' along z direction. A sketch for definition 
of k + and k~ is shown in Fig. 2. The eight coefficients 
a,b, ■ ■ ■ are determined by the boundary conditions (|2j). 
Denoting 

X* = (a,b,A 1 ,B 1 ,c,d,A 2 ,B 2 ), 

with superscript t implying the transpose of vector X, 
we have from Eqs. (O, 

MX = Z, (3) 

where M is a 8 x 8 matrix, and Z is a column vector of 
8 components. Expressing M in terms of 16-block 2x2 
matrices, we get 



M = 



/£>i(0) O(-Z) O(0 

D 2 (<£) rO(l-Z) rO(Z-l) 

rO(l-l) D 2 (-0) rO(l-Z) 

O(0 I>i(-0) O(-Z) 

Dl (<t>) = (: ve t rue t), 

\ rue 4> rve <t> J 



(4) 
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with e± = exp(i/c ± ). The vector Z is given by 

Z l = (-rue^,, -rve^, -ue<pe+, 0, 0, 0). 

Similarly, we can write down the expression for the in- 
coming wave functions ip r from the right superconductor. 
But with the configuration of the SNS junction under 
consideration, ip T can be obtained from the relation 



(5) 



with A = ±1. 

For the wave function ip n of a bound state with energy 
< E n < \A(k±)\, the expression is given by^2 



MJ) = a n (l ne e t)e ik '*'+b n ' U "' 
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FIG. 2: (Color online) Sketch for definition of wave numbers 
k + and k~ on energy curve E(k). The reflected waves of an 
incoming wave of k + include two components of a the Andreev 
and b the normal reflections. 

where k is a complex wave number determined by 
£(fc) = j 7 with 7 = ^/A 2 (kj_)-El (Imfc > 0), w„ = 
exp(i6'/2)/v / 2 with 6* = arctan(7/i£ n ), and <?2 are de- 
termined by £,(q\ l ) + V = -£(q%) -V = E n . The vector 
of the coefficients 

X n = ( a n, b n , A 1 , B 1 , c n , d n , A 2 , B 2 ), 

now satisfies the following equation 

M n X n = 0, (6) 

where M n is a counterpart of M with u, v, fc + , k~, gi and 
g2 replaced with u n , u*, fe, k*, g™ and gj > respectively. 
The energy £"„ is then determined by 

det(M„) = 0. (7) 

The solution to the jth component of X n is given by 
the algebraic complement minor of ljth element of M n 
(multiplied by a factor that is determined by the nor- 
malization condition (VVilV-Vi) = !)■ We note at this 
moment that (a„, b n , A", Bf) = ±(c n ,d n ,A%,B%) at 
(f> = because of ip n (x;(f>) = ±i{j n (—x;—(f>). Therefore, 
for finite phase difference, the coefficients a n , b n , j4", B± 
should have respectively the same orders of magnitudes of 
c n , d n , A 2 , B 2 . There are various approximations based 
on the Andreev and WKB approximation in the existing 
theories! 13 ! 14 ' 15 ! 16 The approximation in Ref. UH corre- 
sponds to b n — c n = A 2 = B" — 0, taking into account 
only the Andreev reflections but neglecting the normal 
reflections at the right NS interface. It is not correct. 
Actually, in a bound state, the electrons and holes are 
bounced back and forth again and again in the normal 
metal. The normal and Andreev reflections at the two 
interfaces are equally important. 

To derive the expression of the current, we start from 
the operator of current density in the continuum model 
of normal metal 

J(x) = -Im[^(x)ViP(x)]. (8) 



For the lattice case, Vip{x) in Eq. {H is replaced with 
[V'O + 1) ~ "00 ~ l)]/2- Taking the statistical average 
with summing up all the contributions from the states of 
positive and negative energies, we obtain 

'-/„ I?'/ 
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(9) 

with 

J f (k) = Re{[A;(0)A_(0) - A* + {-4>)A- (-</»)] sinq, 
-[B* + (<t>)B-(<l>) - B* + (-<f>)B_(-<t>)]smq 2 }. 

where the integral J BZ dk± runs over the first Brillouin 
zone, A± = A\ ± A 2 , and T is the temperature. The 
first term Jf in right hand side of Eq. (JS|) comes from 
the contributions of the free states. The second term is 
due to the bound states. Here, the phase dependence of 
the coefficients of the free waves is explicitly indicated 
by <j> as their argument. Of course, those coefficients of 
the bound states A n, s and B n, s, the energy E n and the 
wave numbers q n, s depend on the phase <f> as well. At 
<t> = corresponding to the equilibrium state, there is 
no current flowing through the junction. The current is 
driven by a finite phase difference. Instead to investigate 
the phase dependence, we here confine ourselves to the 
problem of length L — 21 dependence of the supercurrent 
with fixed phase difference (j> — ir/2. 

We have used two sets of the parameters (t c , A, T)/t 
= (5, 4.75, 0.25) x 1(T 2 (SI) and (6, 3.8, 0.2) x 10~ 2 
(S2) in our calculation. The parameter of electron hop- 
ping through the SN interface was chosen as to = 0.8i c . 
The chemical potential [i and the potential shift Vq in 
the normal metal were set respectively to fi/t = —0.97 
and Vo/t = —0.042, corresponding to the hole densities 
5 S ~ 0.13 in the superconductor and 5 n ~ 0.11 in the 
normal metal (at finite temperature). Shown in Fig. 3 
are the calculated results for the supercurrent as func- 
tion of the distance L = 21 (in unit of c-direction lattice 
constant) between two superconductors. The circles and 
squares correspond to the parameters of (SI) and (S2), 
respectively. For comparison, we also depict the curves of 
J(L) cx exp(— L/£) with £ = 4 and £ = 6. By the proxim- 
ity theory, J(L) should decay exponentially with a much 
shorter coherence length ~ 1 — 2 A for cuprates. Our 
result implies that the supercurrent can flow through the 
junction with L > £ >> £ n . 

To compare the present calculation with the proximity 
theory for the layered system, we here estimate the theo- 
retical coherence length £ c along c axisi According to the 
uncertainty principle, £ c is proportional to the inverse of 
the uncertainty of momentum Sp c of electrons. The lat- 
ter can be estimated as v c Sp c w 8E, where v c is the aver- 
aged magnitude of the electron velocity along c-direction. 
Note that there is no Fermi surface across the c-axis in 
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FIG. 3: (Color online) c-axis supercurrent J through a d- 
wave SNS junction as function of the distance L between two 
superconductors. L is given in unit of c-axis lattice constant. 
The circles and squares are the calculated results for two sets 
of parameters (see the text). The dashed lines express the 
formula exp[— (L — 2)/£] with £ = 4 and £ = 6, respectively. 

the layered system with weak interlayer hopping. From 
the energy dispersion in c-direction, e(q) = — 2t c cos(<7) 
with q the c-axis momentum in unit of c (c-axis lattice 
constant) = 1, we obtain the electron velocity 2t c sin(g>). 



The overall magnitude of v c can be estimated as t c . On 
the other hand, the uncertainty of the energy AE is the 
order of the bandwidth 4t c . We then have £ c « 1/4. This 
£ c is approximately the same as the observed data for 
cuprates. Therefore, according to the proximity theory, 
the supercurrent cannot transfer along the c axis even 
for very short SNS junctions. However, since the super- 
current can be conducted by the Andreev reflections, it 
is not limited by the coherence length. Our calculation 
may account for the giant proximity effect observed in 
the cuprates SNS junctions. 

In summary, we have investigated the supercurrent 
transferring through the c-axis cuprate SNS junctions. 
Due to the Andreev reflections, the supercurrent is con- 
ducted by the in-gap bound states and the free states 
above but close to the gap. It is shown that the super- 
current as a function of thickness of the normal metal 
decreases much slower than the exponential decaying 
expected by the proximity effect. This result implies 
that the supercurrent can transfer through the clean c- 
axis cuprate SNS junctions with the normal metal much 
thicker than its coherence length. The present result may 
account for the giant proximity effect observed in the 
cuprate SNS junctions. 
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